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We present a phenomenologically motivated model which is able to give rise both to an inflationary 
epoch and to the present day cosmic acceleration. We introduce an approach where the energy 
density depends on the scale factor a in such a way that a smooth transition from the inflation to 
the matter dominated era and then to the nowadays accelerated expansion is achieved. We use the 
dimensionless coordinate distance y(z) to the Gold Type la Supernovae dataset and to a sample 
comprising 20 radio galaxies, the shift parameter 1Z and the acoustic peak parameter A to test 
whether the model is in agreement with observations and to constrain its main parameters. As 
an independent cross check, we also compare model predictions with the lookback time to galaxy 
clusters and the age of the universe. Although phenomenologically inspired, the model may be 
theoretically motivated either resorting to a scenario with a scalar and phantom fields (eventually 
interacting) or as an effective description arising from higher order gravity theories. 

PACS numbers: 98.80.-k, 98.80.Es, 97.60.Bw, 98.70.Dk 



I. INTRODUCTION 

An impressive amount of astrophysical data has been 
accumulated in recent years opening the era of the ob- 
servational cosmology and leading to a complete revision 
of what is thought about the universe. The Hubble di- 
agram of Type la Supernovae (hereafter SNela) 0, 0, 
the anisotropy spectrum of the cosmic microwave back- 
ground radiation (hereafter CMBR) HH0, the matter 
power spectrum determined from the clustering proper- 
ties of the large scale distribution of galaxies [6j and the 
data on the Lya emitting regions Q all provide indica- 
tions that the universe have to be described as a spa- 
tially flat manifold where matter and its fluctuations are 
isotropically distributed and represent only about 30% 
of the overall content. In order to fill the gap and drive 
the acceleration, a dominant contribution coming from 
a homogeneously distributed negative pressure fluid has 
been invoked, usually referred to as dark energy. 

Notwithstanding the astonishing flood of papers on 
this topic, it is still unclear what are the nature and the 
nurture of this mysterious component. The simplest ex- 
planation claims for the cosmological constant A thus 
leading to the so called ACDM model 8]. Although be- 
ingthe best fit to most of the available astrophysical data 
the ACDM model is also plagued by many prob- 
lems on different scales. If interpreted as vacuum energy, 
A is up to 120 orders of magnitudes smaller than the 
predicted value. Furthermore, one should also solve the 
coincidence problem, i.e. the nearly equivalence of the 
matter and A contribution to the total energy density. 
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As a response to these problems, much interest has 
been devoted to models with dynamical vacuum energy, 
dubbed quintessence |ll| . These models typically in- 
volve scalar fields with a particular class of potentials, 
allowing the vacuum energy to become dominant only 
recently (see 0] f° r comprehensive reviews). Although 
quintessence by a scalar field is the most studied candi- 
date for dark energy, it generally does not avoid ad hoc 
fine tuning to solve the coincidence problem. 

Relying on the observation that there is actually no 
clear evidence that dark energy is not the same as dark 
matter, many models have been proposed in which a sin- 
gle fluid accounts for both of them. Collectively referred 
to as unified dark energy (UDE) models, these proposals 
resort to standard matter with exotic equation of state 
(hereafter EOS) as the only candidate to both dark mat- 
ter and dark energy thus automatically solving the coin- 
cidence problem. The EOS is then tuned such that the 
fluid behaves as dark matter at high energy density and 
quintessence (or A) at the low ener gy l imit. Interest- 
ing examples are the Chaplygin gas |13j. the tachyonic 
field 0|, the Hobbit model and the Van der Waals 
quintessence |16). 

Actually, there is still a different way to face the prob- 
lem of cosmic acceleration. As stressed in Lue et al. ^3 > 
it is possible that the observed acceleration is not the 
manifestation of another ingredient in the cosmic pie, 
but rather the first signal of a breakdown of our under- 
standing of the laws of gravitation. From this point of 
view, it is thus tempting to modify the Fricdmann equa- 
tions to see whether it is possible to fit the astrophysical 
data with a model comprising only the standard mat- 
ter. Interesting examples of this kind are the Cardas- 
sian expansion 0] and the DGP gravity [l^. Moving 
in this same framework, it is possible to find alternative 
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schemes where a quintessential behavior is obtained as 
a byproduct of fourth order theories of gravity obtained 
by replacing the Ricci scalar curvature R in the gravity 
Lagrangian with a generic function f(R) giving rise to a 
curvature quintessence |20l 12 ll l22l l23l |24| | . 

Despite their radically different underlying philoso- 
phies, the above scenarios share a common feature. They 
all consider cosmic acceleration as a unique event in the 
evolutionary history of the universe and thus focus their 
attention only to a limited redshift range. Actually, 
another period of accelerated expansion took place in 
the very early universe, namely during the inflationary 
epoch. Indeed, it was inflation that has suggested to 
consider a dynamical scalar field as the candidate to ex- 
plain the present day cosmic acceleration. An immediate 
step along this direction consists in wondering whether 
the scalar field today driving the universe is the same 
as the one which originated the inflationary expansion. 
Actually, as scalar fields are not the unique possibility to 
explain the late time acceleration, in the same way they 
are not the only candidates to give rise to the inflation 
with extended theories of gravity as valid alternatives. 

Recently, some attempts have been proposed to explain 
today accelerating universe as a byproduct of primordial 
inflation |25| without invoking any dark energy model, 
but serious objections have been raised [26j making this 
approach controversial. 

From this short overview, it is clear that there is a 
wide confusion among rival theories for both on inflation 
and dark energy. Nevertheless, it is highly desirable to 
find a mechanism to unify both these fascinating issues 
under the same physical scenario. In order to solve this 
problem, we adopt here a phenomenological approach. 
Rather than proposing a somewhat theoretically moti- 
vated model, we introduce a single fluid whose energy 
density dependence on the scale factor is assumed a pri- 
ori in such a way that it gives rise to a smooth transition 
from the inflationary epoch to the matter dominated era 
finally leading to the present day accelerated expansion. 
After having tested the model against observations, we 
also present two possible physical interpreted in terms of 
two radically different scenarios. 

The plan of the paper is as follows. In Sect. II we 
present the main feature of the class of phenomenologi- 
cal models we propose, while Sect. Ill is devoted to the 
study of the dynamics of the universe in presence of this 
fluid. The matching with observations is performed in 
Sect. IV where we consider the dimensionless coordinate 
distance y(z) to the Gold SNela dataset and a sample of 
20 radio galaxies, the shift parameter TZ (related to the 
distance to the last scattering surface) and the acoustic 
peak parameter A. As an independent cross check, we 
also compare the model predictions with the estimated 
age of the universe and the lookback time to galaxy clus- 
ters. Two possible theoretical scenarios to recover our 
model as an effective description are presented in Sect. V, 
while we summarize and conclude in Sect. VI. 



II. THE INFLESSENCE MODEL 

Current theoretical schemes and observational evi- 
dences point out that the evolutionary history of the 
universe comprises two periods of accelerated expansion, 
namely the inflationary epoch and the present day dark 
energy dominated phase. In both cases, the expansion 
is usually interpreted as the result of the presence of a 
negative pressure fluid dominating the energy budget. 
It is natural wondering whether a single (effective) fluid 
may indeed be responsable of both periods of acceler- 
ated expansion. At the same time, this fluid should be 
subdominant during the radiation and matter dominated 
epochs in order to not interfere with baryogenesis and 
structure formation. While it is quite difficult to theo- 
retically formulate the properties of such a fluid, it is, 
on the contrary, clear what are its phenomenological fea- 
tures. Motivated by this consideration, we propose here 
to start from the phenomenology of the fluid in an at- 
tempt to recover its fundamental properties. 

Let us thus assume that the energy density as function 
of the scale factor a reads : 

M-«a-*(l + *)'( 1 + ±)- (!) 

with M a normalization constant, (/?, 7) slope parameters 
and ai << oq two scaling values of the scale factor. For 
later applications, it is convenient to rewrite Eq.Q) in 
terms of the redshift z = I /a — 1 (having set a — I with 
the subscript denoting henceforth quantities evaluated 
at the present day, i.e. z = 0) : 

, (2) ^ (1 + 2)3 ( 1 + i±£) S ( 1 + l^), (2) 

having defined : 

Zj = \ jai - 1 , (3) 

z Q = l/aq - 1 . (4) 
From Eq.Q, it is quite easy to see that : 

{p ~ a~ (/3+3) for a << a/ << a Q 
p ~ a~ 3 for ai << a << ag . (5) 
p ~ a 7-3 for ai << aQ << a 

The energy density of the proposed fluid scales as that 
of dust matter in the range etj << a << aQ. As we 
will see, this means that the fluid follows matter along a 
large part of the universe history, while it scales differ- 
ently only during the very beginning (a << a/) and the 
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present period (a >> a<g). Moreover, choosing (3 = —3, 
the fluid energy density remains constant for a << aj 
thus behaving as the usual cosmological constant A dur- 
ing the early epoch of the universe evolution. Finally, 
the slope parameter 7 determines how the fluid energy 
density scales with a in the present epoch. 

It is still more instructive to look at the equation of 
state w = p/p where p is the fluid pressure. To this aim, 
we first remember the continuity equation : 



3H{p + p) = 



(6) 



with H = a/a the Hubble parameter. Using the obvious 
relation dp/dt — dp /da x da/dt, it is immediate to get 
the following expression for the pressure : 



P : 



i( a t +3p 



(7) 



which holds whatever is the cosmological model provided 
that H is not vanishing everywhere (i.e. the universe is 
not stationary). Actually, if the model contains more 
than a single fluid, this relation strictly holds only for 
the total pressure and the total energy density. However, 
it is still valid for each fluid if each one satisfies the con- 
tinuity equation separately. This happens whenever the 
fluids are not interacting as in many quintessence models. 
It is worth noting, however, that an interaction term be- 
tween the dark matter and the dark energy component is 
an efficient tool to solve the cosmic coincidence problem. 
Nevertheless, as discussed by Zimdhal [27|. the decay rate 
of one substance into the other one enters explicitly only 
in the second derivative of the Hubble parameter and is 
thus quite difficult to be dtected. Moreover, it is also 
possible to show that any cosmic EOS w can be inter- 
preted as an effective EOS of an interacting mixture of 
dust matter and a dark energy component with a suit- 
ably assigned EOS wx- This degeneracy among inter- 
acting and non interacting models ensures that assuming 
that our assumption of non interacting fluids does not 
lead to any loss of generality. Motivated by this consid- 
eration, we may therefore assume that Eq.J7Jl holds for 
each fluid separately and use it to determine the EOS of 
our proposed component. 

Inserting Eq.JTJ into Eq.Q, after some algebra we get : 



P ( l + z 

w = — 

3 \2 + z + Zl 



Z Q 



7 



3 \2 + z + z Q 



(8) 



which is plotted in Fig.^ as function of log (1 + z) for 
models with different values of 7 and zq having set (3 = 
—3 and Zi — 3454 (see later for the motivations of this 
choice). It is worth noting that w does not depend neither 
on 7 nor on zq for high values of z which is expected 
looking at Eq. JSJ . On the contrary, these two parameters 
play a key role in determining the behaviour of the EOS 
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FIG. 1: The equation of state w as function of lz = log (l + z) 
for models with (f3,zi) = (—3,3454). In the left panel, zq — 
1.5 is set and the three curves refer to 7 = 2 (short dashed), 
7 = 3 (solid) and 7 = 4 (long dashed). In the right panel, 
instead, we set 7 = 3 and use three different values of zq, 
namely zq — 0.5, 1.5, 2.5 (short dashed, solid and long dashed 
respectively) . 



over the redshift range (0, 100) which represents most of 
the history of the universe (in terms of time). 

The role of the different quantities (J3, j,zi,zq) is bet- 
ter understood considering the asymptotic limits of the 
EOS. We easily get : 



lim w(z) = — 

:^±00 3 



(9) 



where z —> —00 refers to the asymptotic future. Eq.@ 
shows that setting j3 — —3, the fluid EOS asymptot- 
ically approaches that of the cosmological constant, i.e. 
w\ = — 1. In general, if we impose the constraint (3 < — 1, 
we get a fluid having a negative pressure in the far past 
so that it is able to drive the accelerated expansion oc- 
curring during the inflationary epoch. 

It is now clear that zi controls the transition towards 
the past asymptotic value in the sense that the larger is z 
with respect to zj, the smaller is the difference between 
w(z) and its asymptotic limit /3/3. This consideration 
suggests that Zj has to take quite high values (indeed, 
far greater than 10 3 ) since, for z » zi, the universe is 
in its inflationary phase. 

Let us now consider the present day value of w that 
turns out to be : 



w 



P 7 ( l±*Q 

3(2 + 2/) 3 \2 + z Q 



Z Q 



3 \2 + z Q 



(10) 



where, in the second line, we have used the fact that Zi 
is very large. Being zq > 0, in order to have a present 
day accelerated expansion, wq should be negative so that 
we get the constraint 7 > 0. Moreover, depending on 
the values of 7 and zq , wq could also be smaller than w\ 
so that we may recover phantom -like models 28] . The 
parameter zq then regulates the transition to the dark 
energy - like dominated period. 

Summarizing, the fluid with energy density and EOS 
given by Eqs.Q and © is able to drive the accelerated 
expansion of the universe during both the inflationary 
epoch and the present day period. Therefore, such a fluid 
play the role of both the inflaton and the quintessence 
scalar field so that we dub it inflessence (contracting the 
words inflationary quintessence). 
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It is worth noting that, since p scales with a as the dust 
matter energy density for a long period of the universe 
history, the coincidence problem is partially alleviated. 
Moreover, although wo could be smaller than -1 as for a 
phantom field, the EOS asymptotically tends to a value 
larger than -1 (provided that — 3</3<— 1) so that it 
is likely that the Big Rip is significantly delayed or does 
not occur at all. 



III. THE DYNAMICS OF THE UNIVERSE 

Having worked out the main physical features of the 
model, we are now able to study the dynamics of the uni- 
verse. To this aim, we assume spatial flatness in agree- 
ment with the CMBR observations 3, S ■ As a first 
step, we have to to decide what are the ingredients of 
the cosmic pie. Since we are interested to the full evo- 
lutionary history of the universe, we cannot neglect the 
radiation contribution even if including or neglecting this 
term will have no impact at all on the cosmological tests 
we will perform later. An obvious ingredient is repre- 
sented by the baryons since they make up the structures 
(galaxies and cluster of galaxies and the hot intraclus- 
ter gas) we observe today. It is more subtle to decide 
whether including or not dark matter. As we have seen, 
the inflessence energy density scales as that of dust mat- 
ter along the most of the universe history decoupling from 
it only recently (i.e., for z << zq). From this point of 
view, inflessence should be considered as a candidate for 
a unified dark energy model, similar to the (generalized) 
Chaplygin gas ^3| or the Hobbit model 0] , with the fur- 
ther attracting feature of unifying also inflation with dark 
matter and dark energy. If this interpretation is correct, 
however, the inflessence EOS should be almost vanishing 
over a certain redshift range. Solving w{zm) = 0, we get : 



ZM 



VQl - (2 + Vq)P ± ^y Q Z(l3, 7, z Q , zj) 

2(3 [ ' 



with uq = 1 + zq and 



Z(p, 7) z Q , Zl ) = y Q fi 2 + 2(2 -y Q + 2z I )(3 1 + y Ql 2 . (12) 

It is easy to check that, for reasonable values of /3 and zi, 
Zm is always a complex number and hence the EOS never 
vanishes and is always negative. Thus, even if its energy 
density scales as that of dark matter over the most of 
the universe life, inflessence cannot play the same role of 
matter since its equation of state is always significantly 
different from null. As a result, we have to include also 
the dark matter in the total energy budget. 

Introducing now the critical energy density p crit = 
3Hq /8ttG, the expansion rate H is determined by : 



e 2 {z) = Q r (i + zf + q m (i + zf + n x g{z) (13) 
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FIG. 2: The deceleration parameter q as function of the red- 
shift z for models with (f3,zi) = (—3,3454). In the left pan- 
els, zq = 1.5 is set and the three curves refer to 7 = 2 (short 
dashed), 7 = 3 (solid) and 7 = 4 (long dashed). In the right 
panels, instead, we fix 7 = 3 and consider three different 
values of zq, namely zq = 0.5,1.5,2.5 (short dashed, solid 
and long dashed respectively). Upper panels refer to the case 
fix = 0.7, while the lower ones are for Qx = 0.96. 



with E(z) = H(z)/H , f2j = p t (z = 0)/p cr it the to- 
day density parameter for the i-th component, the sub- 
scripts r, M, X denoting quantities referring to radiation, 
matter (which means the sum of baryons and dark mat- 
ter) and inflessence respectively, and we define g(z) = 
px(z)/px{z = 0). From Eq.lJ^J, we easily get: 

Note that, since we are working in a spatially flat uni- 
verse, we may express Qm in terms of the remaining two 
density parameters. 

Let us now consider the deceleration parameter q = 
-ad/ a 2 = -1 + (1 + z)d\nE(z)/dz. Fig.0 shows q(z) 
for models with (j3, zj) = (—3, 3454) and different values 
of the parameters (7, Zq,Qx)' For this latter quantity, 
in particular, we consider a realistic case (Qx — 0.72) 
and a more extreme one (Qx — 0.96) that allows us 
to investigate the consequences of neglecting the dark 
matter term (so that Qm = Qb — 0.04, being Qi, the 
baryons density parameter). 

As a general remark, we note that the higher is 7 or 
zq, the lower is the deceleration parameter at a given z, 
i.e. the more the universe is accelerating. As a conse- 
quence, the transition redshift zt, defined by the condi- 
tion q(zr) = and marking the transition from decel- 
eration to acceleration, is an increasing function of both 
these parameters. This may be easily understood consid- 
ering that zq controls the transition from the matter - like 
to the dark energy -like scaling of the inflessence energy 
density, while 7 plays the main role in determining the 
value of Wq. Therefore, the higher is zq, the larger is 
the redshift range over which inflessence behaves as mat- 
ter (and hence q > 0) so that zt increases. Similarly, the 
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FIG. 3: The age of the universe to (in Gyr) as function of 
zq for models with (/3,zi,h) = (—3,3454,0.664) and three 
different values of 7, namely 7 = 2 (short dashed), 7 = 3 
(solid) and 7 = 4 (long dashed). Left panel refers to the 
realistic case Qx = 0.7, while the right one to Qx = 0.96. 



higher is 7, the higher (in absolute value) is wo so that the 
universe accelerates for a longer period and zt increases. 
As a cross check, one may study the behaviour of zt with 
the parameters (7, zq) which demands for a numerical 
investigation since the equation q(zx) = may not be 
solved analytically. Actually, such an analysis only par- 
tially confirms what we have said above. Indeed, while 
it is true that Zt increases with zq even if a saturation 
effect appears for zq > 2 3, it turns out that zt may 
also decrease with 7 for values of 7 > 4 -j- 5. However, 
because of the interplay between 7 and zq, it is actually 
quite difficult to infer a general result. 

Not surprisingly, the deceleration parameter is a de- 
creasing function of fix, i-e. the larger is fix, the more 
the universe accelerates. This is quite obvious since the 
higher is fix, the more the negative pressure inflessence 
fluid dominates the energy budget giving rise to acceler- 
ated expansion. As a consequence, zt is an increasing 
function of fix whatever are the values of 7 and zq . It 
is worth noting that, since observations argue in favour 
of zt ~ 0.5 0|, the case fix = 0.96 is clearly disfavoured 
thus suggesting that inflessence may not play the role 
of unified dark energy in agreement with what we have 
inferred above on the basis of the behaviour of the EOS. 

As a final issue, let us determine the evolution of the 
scale factor a with the cosmic time t. To this aim, we 
first integrate H(z) to get the age of the universe at z : 



t(z) = t H 



dz' 



[l + z')E{z>) 



(15) 



where tjj = 9.773/i _1 Gyr is the Hubble time 
with h the present day Hubble constant in units of 
100 km s _1 Mpc -1 . Ea. l|15|) cannot be integrated an- 
alytically in the general case, but it is straightforward to 
compute numerically. In particular, setting z = as the 
lower limit of integration, we get the present day age of 
the universe to that is shown in Fig.|3 From the plot, it 
is clear that to is an increasing function of both 7 and 
zq with the former parameter playing a more significant 
role. On the other hand, to strongly depends also on 
fix with larger values of fix leading to higher values 
of the age of the universe. In particular, this provides 
a further argument against inflessence as a unified dark 



energy model since, for values of fix close to unity, to 
may become unrealistically high. 

Finally, we numerically invert Eq. l|15H and use a = 
(1 + z) _1 to get the dependence of the scale factor on 
the cosmic time. It turns out that the dependence of 
the scale factor on the parameters (7, zq) is more or less 
important depending on the value of fix being quite weak 
for realistic values of fix itself. Only if the inflessence 
plays the role of unified dark energy, the dependence of 
a(t) on (7, zq) is significant. In each case, however, we 
note that, for a given t, the scale factor is larger for lower 
values of 7 and zq. It is worth noting that a(t) does not 
diverge for finite t even if wq is smaller than -1 for some of 
the combinations (7, zq) so that, as yet suggested above, 
the Big Rip is avoided. 



IV. MATCHING WITH OBSERVATIONS 

Whatever model which aim is to describe the evolution 
of the universe must be able to reproduce what is indeed 
observed. It is thus mandatory to test the viability of the 
proposed inflessence scenario by contrasting and compar- 
ing it to the astrophysical data available up to now. This 
is also a powerful tool to constrain the model parameters 
thus paving the way towards a complete characterization 
of the model. As an interesting byproduct, constraining 
the characterizing parameters allows to estimate some 
quantities common to every model (such as go, to, w o, 
zt) to previous values in literature. 



A. The method 

In order to constrain the inflessence model parameters, 
we maximize the following likelihood function : 



L oc exp 



x 2 (p) 



(16) 



where p denotes the set of model parameters and the 
pseudo - \ 2 merit function is defined as : 



A' 



x 2 ( P ) = £ 



y th (z i ,p)-y° i 



TZ{p) - 1.716 



0.062 



A{p) - 0.469 



0.017 



.(17) 



Let us discuss briefly the different terms entering Eq. (|17(l . 
In the first one, we consider the dimensionless coordinate 
distance y to an object at redshift z defined as : 



y{z) 



dz' 
W) 



(18) 
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and related to the usual luminosity distance Dl as Dl = 
(1 + z)y(z). Daly & Djorgovki |2{| have compiled a sam- 
ple comprising data on y(z) for the 157 SNela in the 
Riess et al. Gold dataset and 20 radiogalaxies from 
pjof , summarized in Tables 1 and 2 of [2^| . As a prelim- 
inary step, they have fitted the linear Hubble law to a 
large set of low redshift (z < 0.1) SNela thus obtaining : 

h = 0.664 ± 0.008 . 

We thus set h = 0.664 in order to be consistent with 
their work, but we have checked that varying h in the 
68% CL quoted above does not alter the main results. 
Furthermore, the value we are using is consistent also 
with H = 72 ± 8 km s _1 Mpc -1 given by the HST Key 
project [U based on the local distance ladder and with 
the estimates coming from the time delay in multiply 
imaged quasars (3^ and the Sunyaev- Zel'dovich effect 
in X - ray emitting clusters 33] . 

The second term in Eq. (|17l) makes it possible to extend 
the redshift range over which y{z) is probed resorting to 
the distance to the last scattering surface. Actually, what 
can be determined from the CMBR anisotropy spectrum 
is the so called shift parameter defined as [34], l35| : 



(19) 



where z\ s is the redshift of the last scattering surface 
which can be approximated as |3(j,: 



z is = 1048 (l + 0.00124w, 



-0.738 



)(l + <7lO (20) 



with LOi = Vtih 2 and (51,52) given in Ref. [3y|. The pa- 
rameter ojb is well constrained by the baryogenesis calcu- 
lations contrasted to the observed abundances of primor- 
dial elements. Using this method, Kirkman et al. [37| 
have determined : 

uj b = 0.0214 ±0.0020 . 

Neglecting the small error, we thus set cub = 0.0214 and 
use this value to determine zi s . It is worth noting, how- 
ever, that the exact value of zi s has a negligible impact 
on the results and setting zi s = 1100 does not change 
none of the constraints on the other model parameters. 

Finally, the third term in the definiton of x 2 takes into 
account the recent measurements of the acoustic peak in 
the large scale correlation function at 100 hr 1 Mpc sep- 
aration detected by Eisenstein et al. j3g] using a sample 
of 46748 luminous r ed g alaxies (LRG) selected from the 
SDSS Main Sample 39]. Actually, rather than the posi- 
tion of acoustic peak itself, a closely related quantity is 
better constrained from these data defined as 38] : 



A = 



ZLRG 



ZLRG 



E(zlrg) 



V 2 {zlrg) 



1/3 



with zlrg = 0.35 the effective redshift of the LRG sam- 
ple. As it is clear, the A parameter depends not only on 
the dimensionless coordinate distance (and thus on the 
integrated expansion rate), but also on Qm and E(z) ex- 
plicitly which removes some of the degeneracies intrinsic 
in distance fitting methods. Therefore, it is particularly 
interesting to include A as a further constraint on the 
model parameters using its measured value |38| : 

,4 = 0.469 ±0.017. 

Note that, although similar to the usual reduced \ 2 m ~ 
troduced in statistics, the reduced x 2 (i- e -, the ratio be- 
tween the x 2 an d the number of degrees of freedom) is 
not forced to be 1 for the best fit model because of the 
presence of the priors on 1Z and A and since the uncer- 
tainties o~i are not Gaussian distributed, but take care 
of both statistical errors and systematic uncertainties. 
With the definition l|16l) of the likelihood function, the 
best fit model parameters are those that maximize £(p). 
However, to constrain a given parameter pi, one resorts 
to the marginalized likelihood function defined as : 



dpi-i I dp.. 



Ji+i- 



d Pn C{p) (22) 



(21) 



that is normalized at unity at maximum. Denoting with 
Xo is the value of the \ 2 f° r the best fit model, the 1<7 
confidence regions are determined by Ax 2 = X 2 ~ Xo = 1 > 
while the condition A% 2 = 4 delimited the 2er confidence 
regions. 



B. Results 

Before analyzing the results, it is mandatory to as- 
sess what is the parameter space to explore. In princi- 
ple, when fitting the inflessence model to the observa- 
tional constraints discussed above, we should assign the 
values of eight parameters, namely the slope parameters 
(/3, 7), the scaling redshifts [zi, zq), the present day den- 
sity parameters of the inflessence, radiation and baryons 
(tlx, fir, Qb) and the Hubble constant h. Given such a 
large parameter space, it is quite unlikely that meaningful 
constraints could be derived. Actually, we are interested 
only in a subspace of this 8D space. Indeeds, we have 
already set u>i> = fl^h 2 as explained above, while h has 
been fixed to the same value of Ref. [2{| in order to be 
consistent with their estimate of y (z) for radio galax- 
ies. We also set O r = 9.89 x 10~ 5 |40|. We may further 
reduce the parameter space noting that, since zj is very 
high, both the EOS and the expansion rate (which are 
the quantities probed by the observational data) are in- 
dependent of (5 and Zj. As a consequence, both these 
parameters are unconstrained and may be set a priori 
without affecting anyway the results. Motivated by these 
considerations, we fix f3 = —3 (so that the EOS asymp- 
totically approaches the cosmological constant value) and 
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fix 


wo 


go 


ZT 


to (Gyr) 


bf 


3.73 


0.1 


0.717 


-1.03 


-0.53 


0.65 


14.3 


la 


(3.42, 4.24) 


< 1.9 


(0.698, 0.736) 


(-1.16, -0.89) 


(-0.75, -0.43) 


(0.58, 0.70) 


(13.75, 14.86) 


2(7 


(3.17, 5.86) 


< 5.3 


(0.678, 0.754) 


(-1.52, -0.65) 


(-1.15, -0.19) 


(0.50, 0.75) 


(13.05, 15.52) 



TABLE I: Summary of the results of the likelihood analysis for the infiessence model. The maximum likelihood value (bf) and 
la and 2a range are reported for the model parameters (7, zq, fix) and for the derived quantities (wo, qo, zt, to)- See text for 
details. Note that we are able to give only upper limits on zq whose best fit is at the lower end of the range explored. 
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FIG. 4: Observed and predicted dimensionless coordinate dis- 
tance y(z) diagram for the best fit infiessence model. 



zj = 3454, i.e. the redshift of radiation - matter equiva- 
lence 1 . 

We are thus left with a three parameter space to ex- 
plore defined by the slope 7, the scaling redshift zq and 
the density parameter fix- Applying the procedure de- 
scribed in the previous subsection in this 3D space, we 
finally end with the constraints summarized in Table I, 
while Fig.01 shows the very good agreement between the 
data and the theoretical y(z) for the best fit model. 

It is worth noting that the likelihood analysis makes 
it possible to efficiently constraints both 7 and fix (and 
hence fin), but only upper limits may be given for zq. 
This is also clear looking at the marginalized likelihood 
functions reported in Fig-El This could be qualitatively 
explained considering that zq only controls when the 
dark energy -like behaviour starts dominating. However, 
we only need that Six is larger than flM and the EOS 
is negative today in order to be able to have an acceler- 
ated expansion in the redshift range probed by the SNela 
and radio galaxies. As such, it is not just necessary that 
zq > to fit the data thus explaining why the constraints 



Of course, this does not mean that inflation takes place soon 
before z = 3454. Indeed, the infiessence EOS approaches the 
value w = —1 only for z >> zj. It is thus likely that a more 
correct value should be of order 10 4 ~ 6 , but whatever choice is 
equivalent from the point of view of the effect on the constraints 
on the other parameters. 




Ox 



FIG. 5: Marginalized likelihood functions (normalized to be 
1 at maximum) for the infiessence model parameters 7 (left 
panel), zq (central panel) and fix (right panel). 



on this parameter are so weak. 

In Table I, we also give the estimated values of some 
physically interesting quantities, namely the present day 
values of the EOS Wg and the deceleration parameter 
qo, the transition redshift zt defined by the condition 
q{zx) = and the age of the universe to- Since the 
uncertainties on the model parameters are not Gaussian 
distributed, a naive propagation of the errors is not pos- 
sible. We thus estimate the 1 and 2a constraints on the 
derived quantities by randomly generating 20000 points 
(7, zq,Qx) using the marginalized likelihood functions 
of each parameter and then constructing the likelihood 
function of the derived quantity. Although not statisti- 
cally well motivated, this procedure gives a conservative 
estimate of the uncertainties. 

It is quite interesting to compare our constraints on 
some derived quantities to those previously available in 
literature. First, we consider the total matter density 
parameter that could be derived from our estimate 
of fix thus obtaining fl M = 0.283 ± 0.019 (at la). This 
is in very good agreement with other estimates such as 
fl M = 0.29^03 obtained by Riess et al. by fitting 
the quiessence model (i.e. dark energy with constant 
equation of state) and fl™ = 0.284±°;°™ (at 99% CL) 
found by Seljak et al. [l(| through a combined analysis 
of the CMBR anisotropy spectrum measured by WMAP, 
the power spectrum of SDSS galaxies, the SNela Gold 
dataset, the dependence of the bias on luminosity and the 
Lya power spectrum in the framewrok of the concordance 
ACDM model. It is also consistent with the Allen et 
al. H2 result fl M = 0.245lg;o3° from the fitting to the 
gas mass fraction in galaxy cluysters. Overall, there is a 
very good agreement with our result which is encouraging 
since, whatever is the dark energy component, the matter 
abundance should be the same. 

As a second check, we consider the age of the universe 
to whose best fit value turns out to be 14.3 Gyr in good 
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agreement with to 



in = 14.4 



-1.4 
-1.3 



1 o 24 +0.89 



Gyr from Tegmark et 



Gyr found by Rebolo et al. 



al. 0, 

and to = 13.6±01.19 Gyr given by Seljak et al. 
Although such results are rigorously model dependent, 
the agreement among the different estimates and our own 
leads furhter support to the inflessence scenario. 

It is more interesting to consider the constraints on the 
transition redshift zt- It is difficult to derive this quan- 
tity directly from the data even if some attempts have 
been made. Fitting the phenomcnological parametriza- 
tion q(z) = q + dq/dz\ z= oz to the SNela Hubble dia- 
gram, Riess et al. Q have found z T = 0.46±0.13 (at la 
level) which is in agreement with our result only at the 
2a level. Since it is not clear what is the systematic error 
induced by the linear approximation of q(z), which only 
works over a limited redshift range, we do not consider 
a serious shortcoming of the model the la disagreement 
between this estimate of zt and our one. 

As a final remark, it is intriguing to note that, although 
the 1 and 2a ranges extend in the phantom like region, 
the best fit value of wq is only marginally different from 
that of the cosmological constant. Considering also the 
estimated fix, we may conclude that the best fit inf- 
lessence model is today fully equivalent to the concor- 
dance ACDM model which explains why we are able to 
fit equally well the observational data. 



C. The lookback time to galaxy clusters 

In a recent paper |42l | , some of us have proposed to use 
the lookback time to galaxy clusters and the age of the 
universe as an alternative test to investigate the viability 
of a given cosmological model and to constrain its pa- 
rameters. Being based on a different astrophysics (stel- 
lar evolutionary processes rather than calibrating stan- 
dard candles) and on a different physical quantity (time 
rather than distance), this method is a valid alternative 
to the usual cosmological tests. Since, unfortunately, at 
the moment the dataset is very limited, the method may 
be considered as a cross check of the results obtained by 
the usual probes and in this sense we will apply it to the 
inflessence model. 

Here, we briefly report the basics of the method refer- 
ring the interested reader to Ref. 0] f° r further details. 
First, let us remind that, given an object at redshift z, 
the lookback time is defined as the difference between 
the present day age of the universe and the age of the 
universe at that redshift : 



t L (z,p) = t H 



dz' 



{l + z')E{z') 



(23) 



Having estimated somewhat the age of the object t age , 
t l may be determined from observations as t°^ s (z) = 
to bs _ tage(z) — df , with ig fcs the estimated age of the 
universe and df the so called delay factor which takes into 
account our ignorance of the formation redshift of the 



object (and hence of the difference between the age of the 
universe and that of the object). Constraints on model 
parameters are then obtained by a likelihood analysis 
similar to the one performed in the previous section, but 
redefining the pseudo - \ 2 merit function as 2 : 



Nit 



4-obs 



(*) 



obs 



with <7i and at the uncertainties on the age of the i - th 
cluster and the present day age of the universe respec- 
tively. Following we se t o» = 1 Gyr although this 
is a gross overestimate which does not affect, however, 
the results. We also fix {t° bs ,a t ) = (14.4,1.4) Gyr in 
agreement with |5j. 

Having only seven data (six clusters and the age of the 
universe), it is almost obvious that we are not able to 
constrain the full parameter space (7, zq, fix, df) so that 
we have decided to set fix to the best fit value in Table 
I to finally get the following constraints on (7, zq) : 



3.2 < 7 < 5.4 , zq < 5.1 , 3.1 < df < 4.6 at la 
2.5 < 7 < 7.3 , z Q < 00 , 2.3 < df < 5.3 at 2a 



(25) 



with (7, zq, df) = (4.1, 0.1, 3.8) as best fit parameters and 
df in Gyr. Not surprisingly, the constraints on (7, zq) 
are significantly weakened with respect to the results in 
Table I so that at 2a level the upper limit on zq is for- 
mally unconstrained (i.e., it is larger than the upper end 
of the range probed). However, it is worth noting that 
the constraints on 7 are in very good agreement with the 
ones obtained before and also the best fit value is quite 
close. Moreover, the delay factor turns out to be quite 
similar to that estimated using radically different models 
(see 0] for details) which is a reassuring result since this 
quantity is mainly determined by astrophysical processes 
and should be hence model independent. We may there- 
fore conclude that this test leads further support to the 
result obtained by matching the inflessence model to the 
observational data discussed above. 



V. THEORETICAL MOTIVATIONS 

The key ingredient of the inflessence model we are in- 
vestigating is the scaling of energy density with the scale 
factor given by Eq.Q). This ansatz has been phenomcno- 
logically motivated by the need to have a fluid which is 
able to account for both inflation and the present day ac- 
celerated expansion with the further advantage of track- 



2 Note that in l4Jl there is a prior on h in the definition of \f t that 
we have removed here since h has been set above. 
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ing the matter for a long period thus alleviating the co- 
incidence problem. The very good matching with two in- 
dependent datasets is an a posteriori justification of the 
model and strongly motivates further studies. Neverthe- 
less, Eq.Q) is not based on any underlying physical the- 
ory. Actually given its phenomenological nature, finding 
a unique theoretical background for the inflessence sce- 
nario is not easy. Here, we propose two possible and 
radically different mechanisms which are able to give rise 
to a scenario with a single effective fluid having the same 
properties of the inflessence. 



-0.6 
-0.62 
-0.64 
-0.66 
-0.68 

-0.7 
-0.72 
-0.74 



A. Scalar field quintessence 

Since inflessence is able to fit the observational data 
with the same accuracy as other popular dark energy 
candidates, it is worth wondering whether it is possible 
to recover Eq.([Q) resorting to the same kind of mecha- 
nism. To this end, it is useful to remind that, in most of 
the quintessence models, the cosmic acceleration is driven 
by a scalar field <fi rolling down its self interaction poten- 
tial and giving rise to a fluid whose energy density and 
pressure read : 



p = \i> 2 - W) , 



(26) 



so that the EOS is : 



w<l> = p/p : 



1 - 2V^(0)/0 2 
l + 2V4<t>)/ct> 2 



(27) 



ft is easy to see that < — 1 is not allowed by Eq. (|27ll . 
while this is possible for some combination of the inf- 
lessence model parameters 7 and zq (and, indeed, the 
estimated wq in Table I does not exclude such values). 
Changing the sign of the kinetic term in Eos. (|26|) . we get 
phantom models whose EOS is, on the contrary, always 
smaller than -1. The inflessence EOS may cross the so 
called phantom divide, i.e. the EOS may pass from a 
value w < -1 to w > -1 or vice versa for some combina- 
tion of the model parameters (see, for instance, the long 
dashed curve in the left panel of Fig. [TJ. 

Actually, since zi is very large, for z < zi, Eq.l|5|l may 
be very well approximated as : 



7 

" 3 V 2 



f 



Z Q 



Z Q 



(28) 



and it is easy to show that w > — 1 for z > za = — 1 + 
(f + zq)(1 + 7/3). In order to have acceleration, w; < 1/3 
is needed. Imposing this condition, we get z\ < so that 
the approximated EOS is always larger than -1 and we 



FIG. 6: Reconstructed scalar field potential for the inflessence 
model with parameter values as detailed in the text. Note 
that, for numerical convenience, we plot V Te d = V/p cr it vs 
4> re d = Ho(j>/p cr it so that the numerical values are somewhat 
arbitrary. The reconstruction shown refers to the redshift 
range (0, 3), but nothing prevents to extend it to higher z. 



may recover the inflessence model by resorting to a single 
scalar field. From Eas. l|26|l . we get : 



CO 



V^XPc 

H 



y/[l + w(z')]g(z') 
(l + z')E{z') 



(29) 



v{z) = _^rU {1 _ w{z)]g{z) 



with w(z) and g(z) given by Eqs.© and l|14|) re- 
spectively. Eliminating numerically z from this set 
of equations, we finally reconstruct the potential V((f>) 
which is shown in Fig.H3 for the inflessence model with 
(/3,-f,ZQ,zi,Slx,h) = (-3,3.73,0.1,3454,0.72,0.664). 
As it is clear from the plot, the potential may be very 
well approximated by a simple power law with negative 
slope. Such a shape of the potential is quite common in 
quintessential scenarios. 

For completeness, we remind that, for those combi- 
nations of parameters that give rise to inflessence EOS 
crossing the phantom divide, the interpretation in terms 
of a single scalar field must be abandoned. However, it is 
possible to recover again Eq. (JIJ resorting to the quintom 
scenario where both a scalar and a phantom field are 
introduced with an eventual interaction term. In an al- 
ternative scenario of the same kind, dubbed hessence, the 
role of the two fields is played by the real and the imagi- 
nary part of a single complex scalar field [4^ | . Both quin- 
tom and hessence offer the possibility of recovering the in- 
flessence model for those combination of (/?, 7, zq , zj , Six ) 
leading to w(z) crossing the phantom divide. 



B. Inflessence and f(R) gravity 

Recently, it has been demonstrated that, under 
some quite general conditions, the dynamics (i.e., the 
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FIG. 7: Reconstructed f(R) for the inflessence model with 
parameters set as in Fig.|H| Note that we plot If — log (— /) 
vs IR = log (— R) for numerical convenience and perform the 
reconstruction over the redshift range (0, 3). 

scale factor and the expansion rate) of dark energy model 
may be obtained as the byproduct of a fourth order the- 
ory of gravity. In such a framework, standard matter 
(both dust and radiation) is the only fluid filling the uni- 
verse, but the correct theory of gravity is no more the 
Einsteinian general relativity, but a higher order theory 
obtained by replacing the Ricci scalar R in the gravity 
Lagrangian with an analytic function f{R). It is then 
possible |20l | to show that the dynamical equations may 
still be written in the usual form provided that (i) the 
matter term is coupled to the geometry through f'(R) 
with the prime denoting derivative with respect to R and 
(ii) an effective curvature fluid is defined with energy den- 
sity and EOS determined by f(R) and its derivatives up 
to the third order 3 . 

■ I 

The method developed in [45J makes it possible to de- 
rive the f(R) theory which reproduces a given H (z) and 
may thus be straightforwardly applied to the case of the 
inflessence model. We do not discuss the basics of the 
method referring the interested reader to the quoted pa- 
per, while we show the result in Fig.0 Although an 
analytical solution is not available, we may approximate 
the reconstructed f(R) very well (with an error less than 
l%)as: 

f(R) ~ -dex(/ + filR + f 2 lR 2 + f 3 lR 3 ) (30) 

where dex(x) = 10 x , IR = \og(-R) and (/ , f 1: f 2 , / 3 ) 
depend on the inflessence model parameters. In particu- 
lar, for the choice in figure, we find : 

(/o,/i,/a,/ 3 ) = (-2.15,4.45,-1.92,0.377) . 

We stress, however, that Eq.J3U|) is only a fitting function 
working well over the redshift range considered in the 



3 Because of the presence of this fluid, such scenario is also referred 
to as curvature quintessence. 



reconstruction and it has not to be meant as the actual 
f(R) to be extrapolated to higher z. 



VI. CONCLUSIONS 

The discovery of the present day cosmic acceleration of 
the universe has opened the way to a wide and still on- 
going debate on what is driving such an accelerated ex- 
pansion. Moreover, the evidences in favour of a spatially 
flat and low matter universe coming from the CMBR 
anisotropy spectrum and the large scale structure, re- 
spectively, has made this problem more intriguing. 

Rather than adding a further somewhat theoretically 
motivated proposal to the large list yet available, we have 
here pursued a different approach to the problem. To this 
end, we have first noted that cosmic acceleration is not a 
unique event since the universe has already undergone a 
phase of accelerated expansion in its very early history, 
i.e. during the inflationary epoch. It is clearly desirable 
to have a mechanism able to produce an initial period of 
accelerated expansion followed by a long phase of mat- 
ter domination (not interfering with nucleosynthesis and 
structure formation) finally ending up with the present 
day cosmic acceleration. Motivated by these purely phe- 
nomenological considerations, we have introduced in the 
energy budget a fluid whose energy density is postulated 
to be given by Eq. . By suitably choosing the model pa- 
rameters, the EOS approaches the cosmological constant 
value (wa = —1) both for z — > oo and for z — > — oo so 
that a de Sitter state is recovered in the very early and 
very late universe. Moreover, the fluid energy density 
scales with the scale factor a in the same way as ordi- 
nary (dust) matter for most of the universe history thus 
solving (or, at least, alleviating) the coincidence prob- 
lem. Given its ability to give rise to both inflationary 
and quintessential behaviour, we have dubbed inflessence 
such a fluid. 

A further support to the inflessence scenario is given by 
the matching with observations. Considering the dimen- 
sionless coordinate distance to Gold SNela sample and a 
dataset comprising 20 radio galaxies, the shift parameter 
and the baryonic acoustic peak in the LRG correlation 
function, we have been able not only to show the viabil- 
ity of the model, but also to constrain its main parame- 
ters. The constraints derived on some other interesting 
physical quantities (such as the present day values of the 
EOS, the deceleration parameter, the transition redshift 
and the age of the universe) are in very good agreement 
with previous ones in literature which also comprise some 
model independent estimates. As an independent cross 
check, we have also used a recently proposed method to 
constrain the model parameters with the lookback time 
to galaxy clusters and the age of the universe thus ob- 
taining consistent estimates for the model parameters. 

These results further support the inflessence scenario 
and represent an a posteriori confirmation of the phe- 
nomenological ansatz we are dealing. Nevertheless, we 
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have also wondered whether it is possible to recover 
Eqs.J2J and JSJ by resorting to some theoretically mo- 
tivated mechanism. Two approaches have been investi- 
gated. In the first case, inflessence turns out to be an 
effective fluid resulting from the combination of a scalar 
and a phantom field evolving under their self interac- 
tion potentials with a possible interaction between them. 
However, for the best fit model parameters, a single scalar 
field is necessary and we have numerically reconstructed 
both the evolution of the field and its potential. A rad- 
ically different interpretation of inflessence is connected 
with higher order gravity theories. In such a case, the 
same dynamics generated by the matter and inflessence 
contributions to the energy density is recovered in the 
framework of fourth order theories of gravity. In particu- 
lar, using a recently derived formalism, we have been able 
to numerically determine the function f(R) entering the 
generalized gravity Lagrangian. It is worth noting that 
such an approach allows to find out a particular class of 
extended theories of gravity able to explain both inflation 
and the present day cosmic acceleration. 

The successful results obtained are encouraging and 



motivate further studies of the inflessence scenario. From 
an observational point of view, it is interesting to com- 
pare and contrast the model with the CMBR anisotropy 
and polarization spectrum and with the data on the mat- 
ter power spectrum and the growth index. These tests 
make it possible to check the model over a different red- 
shift range than the SNela and radio galaxies data of- 
fering also the possibility to tighten the ranges for the 
different parameters (in particular zq). It should also be 
desirable to better investigate the theoretical background 
of inflessence. Although two possible mechanisms have 
been proposed, there is still the possibility that other 
approaches may lead to Eq.Q) either resorting to scalar 
fields or some extended theories of gravity and modified 
Friedman equations (e.g., connected with branes). To 
this regard, it is worth stressing that such theoretical 
model should be able to explain at the same time two of 
the most exciting problems of modern cosmology, namely 
inflation and present day cosmic acceleration. As a nat- 
ural consequence, the coincidence problem results to be 
greatly alleviated as it is typical of all unified models. 
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